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Abstract 

the generalized Einstein action is treated quantum 
mechanically by using a quadratic lagrangian form. The canonical 
quantization of this action is obtained by using the auxiliary variable 
to define the generalized momentum. Physical constraints are imposed 
on the surface term, which is defined to be the cosmological constant. 
One obtains the familiar Wheeler-de Witt equation. The solution of 
this equation is in conformity with General Relativity (GR) .In addition 
to the fact that it is free from GR setbacks at the early universe, since it 
gives time decaying cosmological constant. 

The wave function of the universe and the cosmic scale factor 
are complex quantities, which indicates the existance of quantum 
effects. 2. 7k ,CBR temperature is calculated . 



l.Hamiltonian of curvature squared action: 

Palatini formalism which is based on independent variations of the metric and 
the connection, has been known to be equivalent to the corresponding metric 
formulation for the lagrangian which depends linearly on the scalar curvature 
constructed out of the metric and the Ricci tensor of the connection .In this case, in 
fact, field equations imply that the connection has to be the Levi-Civita connection of 
the metric, which in turn satisfies Einstein's equations [1]. This means that there is a 
full dynamical equivalence with the corresponding metric lagrangian , which a 
posteriori coincides with the Hilbert lagrangian of general relativity . 
A generic fourth -order theory in four dimensions can be described by the action [2] 

A= j f(R)^gd n x 1 - 1 

M 

Where /(i?)is a function of Ricci scalar R. 
The field equations are: 

f{R)R aP - \f{R)g a p = f(Ry ap {g ai i 9p v - g a p g^) 1-2 

This can be expressed in the more expressive form 

G a p = Rap ~ 2 9a P ^ = ^ap 1 — 3 

T a p = y a p (/00 - Rf(.R) + f(.Ry aP gp v - g a p )] 1 - 4 

Where T a p is the curvature-energy tensor. The prime indicates the derivative with 
respect to R. 

It is possible to reduce the action to a point-like, Friedmann-Robertson-walker(F R 
W) one. We have to write 
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A = j dtL(a,a,R,R) 



Where dot indicates derivative with respect to the cosmic time. The scale factor (a) 
and the Ricci scalar R are the canonical variables. This position could seem arbitrary 
since R depends on a, a, a. The definition of R in terms of a, d, a introduces a 
constraint which eliminates second and higher order derivatives in action(l-5),and 
gives a system of second order differential equations in {R,a}. 
Action (1-5) can be written as: 

(ad 2 
a a 1 a / y 

Where the lagrangian multiplier A is derived by varying action (1-6) with respect to R 
.this yield 

A = a 3 f(R) 

If we consider homogeneous and isotropic Robertson- Walker metric in form of 

dr 2 

+ r 2 (d9 2 +sin 2 6d<p 2 ) 



-J 



dt\a 3 f{R)-A 
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dS 2 = g^ v dX^dX v = dt 2 - a 2 (t) 
For which 
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I a a k > 

Rjj = I — I T H T 



9ij 



\a a*- as- f 

R = g ij Rij 

Then scalar curvature R will be written as 
( d d 2 k\ 



R = 6\ 



a a 2 a 2 



If we assume that f(R) = Ri, action (1-6) will be written as: 

A = 2(1 -n)n 2 p j a 3 
If we assume / = 4, 



l R2dt 



, f „ fa a 2 k\ 
A = -216n 2 B a 3 - + + - T ) dt 
J \a a L a 1 ) 
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dt 



1-12 



If we put m = —216n 2 p 

ff k 2 d 4 d 2 ^ 

A = ml [ad 2 + 2d(d 2 + k) + — + — + 2k — 
J \ a a a j 

The corresponding Lagrangian in a minisuperspace model involves field 
variables along with their first derivatives and second derivatives. Therefore, in order 
to cast the action in canonical form that demands the corresponding Lagrangian be 
expressed in terms of the field variables along with their first derivatives only, one 
requires to introduce auxiliary variables[3,4,5,6,7,8,9]. It has been observed that any 
form of auxiliary variable produces correct classical field equations in homogeneous 
and isotropic model. Further, it is not enough to get the correct classical field 
equations only, rather one should also be able to produce correct and well behaved, 
quantum description of the system under consideration. By the above statement we 
mean to find a variable, if it at all exists, with respect to which the Hamiltonian would 
turn out to be hermitian. 

Starobinsky [10], was The first one who established the relevance of the fourth 
order gravity theory in cosmology though it appeared earlier in the context of 



2 



quantum field theory in curved space time. Starobinsky [10] presented a solution of 
the inflationary scenario without invoking phase transition in the very early Universe, 
from a field equation containing only the geometric terms. However, the field 
equations could not be obtained from the action principle, as the terms in the field 
equations are generated from the perturbative quantum field theory. Later, 
Starobinsky and Schmidt [11] have shown that the inflationary phase is an attractor in 
the neighborhood of the fourth order gravity theory. 

Quantum cosmology is most elusive as the role of time is not unique and one 
is being unable to define the Hilbert space[12]. This is due to the fact that 'time' is not 
an external parameter in general theory of relativity, rather it is intrinsically contained 
in the theory, unlike its role in quantum mechanics or quantum field theory in flat 
space time. In curved space time, different slices correspond to different choices of 
time leading to inequivalent quantum theories. Likewise, the canonical quantization of 
gravity is devoid of an unique time variable and hence the definition of probability of 
emergence of a particular Universe out of an ensemble is ambiguous. The canonical 
quantization of Einstein-Hilbert action together with some matter fields yields the 
Wheeler-deWitt equation which does not contain time a priori, although it emerges 
intrinsically through the scale factor of the Universe. However, if the canonical 
variable is so chosen that one of the true degrees of freedom is disentangled from the 
kinetic part of the canonical variables, then this kinetic part in the corresponding 
quantum theory yields a quantum mechanical flavor of time. This is possible only if 
the Einstein-Hilbert action is replaced by curvature squared action or modified by the 
introduction of curvature squared term, in the Robertson- Walker minisuperspace 
model. 

Such wonderful relevance of higher order gravity theory in the context of cosmology 
inspired some authors to give interpretation of the quantum cosmological 
wavefunction with R 2 [13] and even R 3 [14] terms in the Einstein-Hilbert action. 
Further the functional integral for the wavefunction of the Universe proposed by 
Hartle-Hawking [15] runs into serious problem, since the wavefunction diverges 
badly. There are some prescriptions [16] to avoid such divergences, though a 
completely satisfactory result has not yet been obtained. However, to get a convergent 
functional integral,Horowitz [9] proposed an action in the form 

S = j d 4 Xj7j [ACf jkl +B(R- AX) 2 ] 1-13 

Where Cfj kl is the Weyl tensor, R is the Ricci scalar, X is the cosmological constant 
and A,B are the coupling constants. 

The action (1-13 ) reduces to the Einstein-Hilbert action at the weak energy limit. To 
obtain a workable and simplified form of the field equations, one may consider a 
spatially homogeneous and isotropic minisuperspace background, for which the Weyl 
tensor trivially vanishes. A.K.Sanyal and B.Modak [3], considered the action like (1- 
11) retaining only the curvature squared term. The field equations for such an action 
can be obtained by the standard variation principle. 

In the variational principle, the total derivative terms in the action are extracted and 
one gets a surface integral which is assumed to vanish at the boundary or the action is 
chosen in such a way that those surface integral terms have no contribution. However, 
for canonical quantization this principle is not of much help and one has to express 
the action in the canonical form, which is achieved only through the introduction of 
the auxiliary variable. Auxiliary variable can be chosen in an adhoc manner, and 
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different choice of such variable would lead to different description of quantum 
dynamics, keeping the classical field equations unchanged. In view of this 
Ostrogradski [7] in one hand and Boulware etal [6] on the other, have made definite 
prescriptions to choose such variables. Ostrogradski' s prescription [7] was followed 
by Schmidt [17]. 

in Boulware etal's [6] prescription the auxiliary variable should be chosen by taking 
the first derivative of the action with respect to the highest derivative of the field 
variable present in the action. 

Hawking and Luttrell [18] utilized Boulware etal's [6] technique to identify the new 
variable and showed that the Einstein-Hilbert action along with a curvature squared 
term reduces to the Einstein-Hilbert action coupled to a massive scalar field, assuming 
the conformal factor. Horowitz on the other hand, showed that the canonical 
quantization of the curvature squared action yields an equation which is similar to the 
Schrodinger equation[9]. Pollock [19] also used the same technique to the induced 
theory of gravity and obtained the same type of result as that obtained by Horowitz 
[9], in the sense that the corresponding Wheeler-de Witt equation looks similar to the 
Schrodinger equation. 

The striking feature of Boulware etal's [6] prescription is that it can even be applied 
in situations where introduction of the auxiliary variable is not at all required, e.g. in 
the induced theory of gravity, vacuum Einstein-Hilbert action etc. 
The classical field equations remain unchanged with or without the introduction of 
the auxiliary variables. 

2 Hamiltonian with respect to ( a): 

Define the auxiliary variable according to Horowitz [9] 

nab ab ^ 9 — 1 

Where h ab is the three metric. 

Action(l-5) can be written in the following form: 

k 2 d 4 d 2 \ 
ad 2 + 2d(d 2 + k) + — + — +2k — )dt 2-2 
a a a J 

Where the Lagrange is 

L = m [ad 2 + 2d(d 2 + k) + — + — + 2k — ) 2-3 
y a a a J 

The auxiliary variable Q is defined as the variation of the action with respect to the 

higher derivative term d ; 

dA (d d 2 k\ n 

— = 2m - + ^ + ^ U 2 2-4 

da \a a / a 1 J 

According to [3], every action should be supplemented by an appropriate boundary 
term. Here we rid of all the zero-order derivative terms to the introduction of auxiliary 
variable. 

k 2 

The zero-order derivative term in action (1-16 ) is — . 

a 

Then action (1-16 ) is written in the following form: 

f ( d 4 d 2 \ f k 2 

A = ml [ad 2 + 2d(d 2 + k) + — + 2k — )dt + m \ — dt 2-5 
J \ a a J J a 

Now we can write action ( 2-2) in canonical form as: 



A = m j 
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f ( O 2 k 2 \ f k 2 

A= \ \—^—-m — )dt+m\ —dt 2-6 
J \4ma a J J a 

f a 2 +k 

y = -a- dt 2-8 

J a 

rk< 

a — m — 
J a 



k 2 

dt 2-9 



ydt 

Finally, the action became in canonical form with the suitable boundary term like 

W(-^ + =r"4) £it+El+ ' 7 2-10 

Assuming all boundary terms vanish at boundary E x + a = 0. This leads to: 
k = 0,and the scale factor (a)is a constant. 

Now we apply Eular-lagrange equation (^|^~|^-=o ) to the action to recover definition 
of the auxiliary variable 

fad 2 k\ 

Qi = 2maa + 2m(rir + k) = 2m - + ^r + ^- cr 2 — 11 

\a a 1 a 1 J 

This is the same definition of the auxiliary as it was defined before. 
Now it will be easy to obtain the Hamiltonian according to: 

H{q,p,t) = diPi -L{q,q,t) 
dL 

dL 



Q 2 mk 2 

H = YPy — + 2-12 

For canonical quantization we express the Hamiltonian in terms of the basic variables, 
for which we choose 

dL 

P X =Qi=~ 
x = d 
PQi = -Y 

then we obtain the Hamiltonian in terms of the basic variables as: 
p 2 mk 2 

H = fPy--p^ + 2-13 

1 4ma a 

As we know H = then; 

p 2 mk 2 

^ Y 4ma a 
To get the dynamical equation( p = —ihV ) 

a h 2 d 2 mk 2 

ih lT = -A "TT ~\ ~ 2-14 

°y \may ax L ay 

This can be expressed as: ih^- = 

Where: H is the effective Hamiltonian. 



_ h 2 



Amay dx 2 
mk 2 



+ V el 



el 



ay 



2-15 
2-16 



Represents the effective potential. 

Assuming perfect fluid, which obey the equation of state p = 00 p 
now we can solve for p and p in equation (3-9) and with aid of equation 

h a + 2 &+^ = 4nG(p-p). 

If the potential (2-16) remains constant with respect to time (t).we found the relation 
of energy density and the scale factor , 



d 

dt\ a r 



= 



ay + ay = 

, f a 2 +k _ 

d z + d dt + ad + or + k = 

J a 

Add k and divide by a 2 , after all we find: 

2d 2 2k a a f a 2 + k k 
— 5- H — 7 H — = r dt H — 7 

a L a 1 a a L J a a L 

Now we can substitute equation (2-17) into {± + 2g) 2 + 2^ = 47rG(P - p)j 

a f a 2 + k k 

AnGp{l -oi) = -—\ dt + -r 

a L ) a a L 
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P = 



at 



a' 



47rGp(l — a)) = 5- 



a 3 



W / J a a 



P = _ 
P 

CO! = 



d + (^-2)J^dt + ^ 



a a 
d 2 +k 



a + r™- 2 )f— dt + ™ 



- 1 



2-18 



2-19 



2-20 



2-21 



we consider two cases as: 

1. a -> 00; d « l;k * 
This 

p d 
-~ - 3- 
p a 

or 



implies 



that 



2-22 

2-23 
2-24 

2-25 



p oc a~" 
p~0 
with 

This represent matter dominated universe, this means that our model described matter 
era at large radius and slow expansion rate, also we found that the density is 
decreasing with increasing of the radius in conformity with GR and experiments. 
2. a -» 0; d -» 00; k 
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p a 

-~-2- 2-26 
p a 

where 

\ 2-27 



which is requires negative pressure 

p~ - \p 2-28 
This result represents early universe ,which predict inflation. 



3. Hamiltonian with respect to ( a 2 ): 

Now according to Boulware's prescription[6], Hamiltonian formulation of an 
action containing higher order curvature invariant term requires introduction of 
auxiliary variable, 

And the extrinsic curvature k ab is: 

hi = k 2 2 = ^33 = -o-a. 3-1 
The action of R 2 can be written as: 

f ( k 2 a 4 d 2 \ 

A = m\ [ad 2 +2d{d 2 +k) + — + — +2k—\dt 3-2 
J \ a a a J 

k 2 

If we get rid the zero-order derivative — to auxiliary variable [4], 

v4 A 2\ r 1,2 



[( 9 ? o 4 d 2 \ fk 2 

A = m I I ad 2 + 2a(d 2 + k) + — + 2k— 1 dt + m I — 

A=m S{ 



dt 3-3 



t „ a. a 2 \ 

ad 2 + 2a(d 2 + k) + — + 2k — )dt + a 3-4 
a a I 



Now we supplement the action with boundary term 

a = m — dt 



= m — i 
J a 



This action can be written in the following canonical form, with 

dA fa d 2 k\ 
Qi = tt = 2m - + —r + ^r)a 3-5 

a 2 \a a 2 a 2 



A = j {-Q 2 d 



aQ 2 2 k 2 \ 

+ m—)dt+l 2 +a 3-6 

4m a 



a = —aa 



j(k)dt 



Assuming all boundary terms vanishing at boundary 1 2 + a = . This leads to the 
same result as before: 

k = 0,and the scale factor (a)is a constant. 
As before we recovered the definition of the auxiliary variable by applying Eular- 
lagrange. 
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For canonical quantization, we express the Hamiltonian in terms of the basic 
variables, for which we choose: 



t-J* 

' • 2 

x = cr 



dt 



and 

Qz = Vx 

Now we obtain the Hamiltonian(H = 0) as: 



-aPa = 



aQ 2 



+ 



mk' 



4m a 
We express the dynamical equation as: 

^ d ah 2 d 2 ^ mk 2 

da Amd dx 2 ad 
_ ah 2 d 2 



where: 



\md dx 2 



3-8 

3-9 
3-10 



V e2 = 



mk L 



aa 



3-11 



4nGp(l -a)) = - _ 3 ) | kdt + 3k j 

47iGpCl-co) = {^}[-jkdt+kl] 

p_- _4@- 3 ) /fedt+3fe jl 



-5 kdt +kt 



0) 2 = 



(Jf - 3) / kdt + 3/c| 



- 1 



3-12 
3-13 
3-14 

3-15 



According to re 
equivalent : 

1. a -> 00; d « l;fc * 



this 
<o 2 ~0 
and 



- f kdt+ki 
J a 

evant treatment of the potential, also, there are tow cases, but they are 



p a 
-~ - 3- 
p a 



implies 



that 



3-16 



2. a 0; d 00; /c gt 



p~0 
p oc a 



-3 



this 
<o 2 ~0 
and 



p a 
-~ - 3- 
p a 

implies 

p~0 
p oc a" 3 



that 



3-17 
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4. Comparison with A.K Sanyal work: 

a 2 

A.K Sanyal (according to Boulware's prescription) choose a variable z = — then he 

found the definition of the auxiliary variable Q as: 

dA _ fa a 2 k\ 

Q = - =liin 2 pi _ + _ + _y 

Then Sanyal suggestion is to get rid of all the total derivative terms prior to the 
introduction of auxiliary variable. After doing this, he found the boundary term that 
supplements the action is different from that one which supplements the action in 
references [3,4,5] when he followed Horowitz[9].A particular action must be 
supplemented by the same boundary term, independent of the choice of variable, 
we found the same boundary terms in this work according to both Horowitz and 

c k Z 

Boulware's prescription, o = m J — dt 

At boundaries we found k = 0,and the scale factor (a)is a constant. 



5. Acceleration and deceleration scenario: 

In this model we did not include matter term in action (1-12), so, this model 
represents curvature effects in all universal parameters, which will be added to that 
parameters which was obtained from matter lagrangian. 
The total effective Hamiltonian of/? 2 in our minisuperspace model is : 

^ ^ ^ h 2 d 2 nh 2 d 2 

H = H +H =JL_" + ^" +V 5_i 

Amav ox' Ama ox' 



mk 2 rl 1 

v = v el+ v e2 =— ? + 5 
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a 

The wave function is: 

2imkx -2imkx 2ikm x —likmx 

xp = xp 1 + ip 2 = C 1 e~lT +C 2 e~ ~ + D 1 e~^~ + D 2 e ah 5-3 

co = ^ " a -^ + V3fl2 - JJ 5-4a 

_ 3 J^ dt + 3 ^ -j kdt + k a 

J a a " 

d-(q + 2)5 + 3- ( q -+l)t-H- 1 
a) = r ^ + ^ L —, 2 5-46 

a 

Where q = — — ,the deceleration parameter, H is Hubble parameter H = - ,B = 

r—dt 

J a 

There are two approaches, first of them represents the early universe: 
1. a^O; 

at this early time we found that, a>, is dominated by the term: d — 2 J — dt 

,to study how this term leads to acceleration universe, write friedmann 
equation: 

d AtzG 

-= =- p(l + 3o)) 

a 5 

c a 2 

• If d « J —dt ,this case happen when we assume that the universe has phase 
transition, which change d from zero to infinity, we can rewrite the integral; 



[a 2 [a 
— dt = - 
J a J a 



da 



now we find the upper critical value of co: 



co —* — , and, a -> .this means that, curvature did not involve in universe 
revolution. 

• If d is less than but comparable with J — dt , we find: 

co < — - ,and, a > 0. this implies negative pressure,p < ,and p > a 
The negative pressure leads to inflation [20]. 

Now if we let a oc t n ,we find the following relation: 

3 - 4n 

co = n 5 — 5 

3n 



1Q n 

e- 

A: 

2- 



0" 
-2- 

-6; 

-10 : 



0.2 



0.4 



0.6 



0.8 



n 



fig(5 — l)variation of oo with respect to n, a -» 



From above equation when 
• |n|»l,ft)->-|. 



• Inl « 1, <u -> - 

n 



then, 



4 1 

--< to < — 

3 n 
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2. Second approach :a -» oo; d « 1; 

• ~| ~ ,this yields one value of a>, where a> = , and, d < ,p oc a -3 ,p = 

0,this represents matter(relativistic particles) , 
If we turn back to action (1-6), adding matter lagrangian, means, adding a pressure 
term [ 21]. So, according to second approach the universe will contract very fast. May 
be back to the beginning and, pre-big bang (oscillating). 

6. Cosmological constant: 

The boundary terms in equation (2-10) and (3-6) represent the cosmological constant. 

k 2 

A 1 = Q 1 y + m — 6 — la 

a 

k 2 ( [a 2 +k \ 

Ai — m 2m a + \ dt (3dd + ad) 6 — lb 

a \ J a J 

Q 1 = 2m(3dd + ad) 
If the cosmological constant equal zero: 

A 1 = 
k 2 

Q 1 Y + m— = 
a 

For simplicity we let 

k = 

In this case 

I — d — I -dt (3dd + ad) = 6 — 2 

Solving equation (6-2) 

3dd' + ad = 6 — 3 

f d 2 

d + — dt = 6-4 
J a 

integrate equation (6-3) : 

ad + d 2 = C 

d 2 a /da\ 2 

a— ^+ — + C = 6-5 
dt 2 \dt) 

Differentiate equation (6-4): 

d ( fa 2 \ 
Tt{ d + j a dt j = ° 
ad + a 2 = 

d 2 a /da\ 2 

a 5? + y =° 6 - 6 

If the integration constant C = 0,equation (6-5) construe to equation (6-6). 

Solving equation (6-6),if a oc t n or a = At 71 , A is a constant. 
We find that: 

t 2n " 2 [(n- 1) +n] = 
1 

n = 2 

By substituting in equation( 4.3-lb) that a = At 71 we find: 

f d2+k dt - An2tn ~ 1 _ fetl ~ n 

■* a n-1 ^l(n-l) 
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Ai = 



mk 



2m 



nA 2 t n (n-l)(n-2)t n - 3 ) 
Now from equation(3.3-6) 

k 2 

A 2 = Q 2 d + m — 



Ant 71 - 1 - + ^—] (3A 2 n 2 t n -\n - l)t n ~ 2 + 

A(n-1) n-1 J 

6-7 



a 



A? — m 2m 

a 



a + 



2aaa — ak — cr 



(ad + j kdtj 



6 - 8a 
6 -8b 



aa 



a + 



A 2 = 
k = 
2add — a 3 



a L 



= 



ad — 

2add — a 3 
a + = = 



6-9 
6-10 

6-11 



d H = C 

a 

When C = 0, equation (6-11) leads to the same differential equation as equation (6-6). 
By substituting in equation(6-8b) that a = At 71 we find: 

TTtk ^ 

-2m An(n- l)(n-2)t n " 3 



A, = 



At 11 



2A 3 n 2 t n t n - 1 (n - l)t n ~ 2 - Ant n ~ x k - A 3 n 3 t 3n ~ 3 



A 2 t 2n 



{A 2 nt n ~ 1 t n 
6-12 



Hot 



A 1 +A 2 



+ fct) 

The total cosmological constant is 
2m 

A tot = — {-10 A 4 n 4 t 3n ~ 4 + 13A 4 n 3 t 3n ~ 4 - 2kA 2 n 3 t n ~ 2 + 10kA 2 n 2 t n - 2 

-4A 4 n 2 t 3n ~ 4 + nk 2 t~ n - 4kA 2 nt n ~ 2 + k 2 t~ n } 6-13 

1 

A = when n = or — with k = 
Cosmological constant vanishes at radiation era. 

Now we can calculate the density parameter 

A _ At 2 

° A " 3H 2 ~ 3n^ 

2mt 2 

H A = — —r{-10 A 4 n 4 t 3n ~ 4 + 13A 4 n 3 t 3n - 4 - 2kA 2 n 3 t n ~ 2 + \0kA 2 n 2 t n ~ 2 
3An L 

- AA 4 n 2 t 3n ~ 4 + nk 2 t~ n - AkA 2 nt n ~ 2 + k 2 t~ n } 6-14 

m 

p A = -^^{-ioA 4 n 4 t 3n - 4 + Y3A 4 n 3 t 3n ~ 4 - 2kA 2 n 3 t n ~ 2 + 10kA 2 n 2 t n - 2 
AAnG 

- AA 4 n 2 t 3n ~ 4 + nk 2 t~ n - 4kA 2 nt n ~ 2 + k 2 t~ n } 6-15 
Now if we study the behavior of A with respect to n ,k and t we find: 



1. n^O, 



A n =o = 



2mk 2 



Pa = 



mk 2 
An AG 



6-16 
6-17 



12 



n A = 



2mk' 



t 2 = 



2mk' 



6-18 



A 3An 2 ' 3A 
From above, when n -> 0, cosmological constant and vacuum energy density, remain 
constants. But the density parameter, vary with time as: 

• t -> oo, n A -* cxi. this seems to be rejected, because n tot = 1. 

• t ~ n, n A 



2mk' 



3A 



A 







1 



n = - ,k 

2 



L n=0.5 



3mk 2 mkA 



20- 



10- 




n=1/2,k=-1 

fig(6 — 1): shows comparison between cosmological constants , n = - and n 



= ~,(/c = -D, 



13 



0.2 



0.4 



0.6 



0.8 



n=2/3 
n=1/2 



•f 

fig(6 — 2): shows comparison between cosmological constants , n = - and n 

3 4 

4mk 2 t2 2mkAti 
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7. wave function and quantization: 

Now we can construct the wave function from the above Schrodinger-like equation. 
Equation (2-15) is second order differential equation (H=0): 
mk 2 

7-1 



h 2 d 2 ip . , ., 

a HTT H —4> = 



2imkx 



-2imkx 



'ip 1 —Q. 1 e h +C 2 e h 

where 

x = d 

The constructed wave function from second Hamiltonian(3-10) 
mk 2 



7-2 



ah 2 d 2 ip 
Amadx 2 



+ 



aa 



ip = 



is: 



where 



2ikm x —likmx 



7-3 



7-4 



dda 



Cx.C-i, D x and D 2 are constants. 
At boundaries, (a) is constant , xpx = ip 2 = 0. 
Since the wave function depends on k(k = 1, —1,0) and 
d(d > 0). 

So k controls the wave function, as the following: 

1. k = 0: 
xp 01 =C 1 +C 2 
and 

xp 02 =D 1 +D 2 

2. k = l: 

-2imkx 

ip +1 =C 2 e~h~ 
and 

—2ikm x 

ip +2 = D 2 e ah 

3. k = -1: 

2imkx 



2 ikm x 

xp_ 2 = D x e ah 
Now if we assume that a = At n 

2imkx / 2n \ikmx -2imkx 

ip(x) = CiQ~~h~ + D 1 ey2^)~^~ + C 2 Q~ ir ~ + D 2 e~y^i 

k = l 

-2imkx ( 2n \ ikm x 



7-5 
7-6 

7-7 

7-8 

7-9 
7-10 



2n \ikm x 



; \2n-l) 



7-11 



/ 2n \ikm 
+ D 7 e~\2n-\)~h 



ljj + =C 2 e h T U 2 

k = -l 

2imkx / 2n \ikm x 

Xp_=C l Q~ ~ + D 1 e^2n-l)^~ 

From boundary conditions when: 

• x = (this is equivalents to n = 0, we found that: 

C 1 +D 1 + C 2 +D 2 =0 7-12 

• n = - and k = 1 ,we find that D x = O.if k = — 1 we find D 2 = 0. 
So, the general form of the wave function is: 
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2i|fc|mx -2i|fc|mx ( 2n \i\k\mx 

ifj(x) = Qe - ft - + C 2 Q~^~ - (C x + C 2 )e"toJ^^ 7-13 
We can recognize that from(7 — 5,7 — 6,7 — 12) flat universe,/c = 0, is equivalent to 
static one a = constant . 

1. n = 1, a oc t ,we find the second wave function phase reduced to the first one 

2i|fc|mx -2i|fc|mx 2\k\m.X 

xp(x) = Cie - S~ -de - S - = 2idsin- L ^ — 7-14 

i i 

2. As we know at radiation era a oc t2 , n = - ,then, the wave function takes the 
following form: 

2i|fc|mx -2i|fc|mx 

xp(x)=C l Q~^~ + C 2 Q~^~ 7-15 
We get from 7-14 and 7-15 C x = -C 2 .D 2 = 

Using 7-14 to find out d by normalizing the wave function, j xpxp*dx = 1 



/ 

o 



A 2 



2\k\mx 
2i Ci sin — - — 



dx = 1 



d = | 7-16 

1 -|sin 2Ab-2Ab 

Where . b = i^L . 

ft 

According to this limit — oo < - < — 1 and since A > , Ci is complex 

number . 

We can write 7-14 as: 



xp(a) = 2 



sinM 7 — 17 



^ 2i46 - sin 2Ab 

Also, we found that 

< n < 1 7-18 
According to 5-5 and 7-18 , at early universe vacuum dose not dominate the universe. 
So the limit of a> is reduced to: 

1 1 

-- < a) < - 7-19 
3 n 

We can find the probability, which we can call it density parameter fi, to our universe, 
from 7-17 

P = n = — , — sin 2 (6d) 7 -20a 

2Ab - sin 2Ab K J 

8\k\m n (2\k\m \ 

P = a = — ir^sin 2 -V— d 7-206 

4A\k\m-hsm 4 -^ V h J 

ft 

When 2Ab » sin 2Ab , we find that H = -sin 2 (M). 

A 

If we take the probability in flat and unaccelerated universe(appindixI-3), 
3 

lim n = - 7-21 

k=0,d=A A 

Now ,if we take the limit of the probability in flat and accelerated universe, 
3d 2 

limn = ^ 7-22 

fc=o A 3 

Using Friedmann equation ,where k = we find: 



16 



limn = 

fc=0 



8nGpa" 
A 3 



AAm — ft sin 
8m 



AAm. 



Maximum probability, £1 = 1 

■ h ■ -i 
a c = - — sin 

2m 

Where a oc t n . 

ft 



occurs 



when, 



a c = A 



2Anm 
ft 



■sin 



-l 



AAm — ft sin 



8m 



n-1 



2Anm 



sin 



AAm— ft sin 



4/4m 



8m 



n-1 



lim a c = A 

n->0 

lim a c = oo 

n-»l 

d r 



2i4m 



■sin 



-l 



AAm— ft sin 



8m 



1-n 



7-23 
* 0, 

7-24 



7-25 



7-26 



7-27 

7-28 
7-29 

7-30 



and: 



lim H r = — 

The subscript c) refers to critical value ,which occurs at maximum probability. 
From 7-20, fi = -jjpr , we find: 

3|/c|m H 2 _ /2|/c|md\ 



P = 



7rG 



[4i4|fc|m- 



ftsin 



4,4|fc|m~ 



sin' 



Where k ^ 0, we find: 

3m tf 2 

P = 



nG 



\4Am — ft sin 



nvn-i/ 



Pc = 



8?rG 



AAm' 



ft 



sin 



/2md\ 
\~ft~J 



7 - 31a 



7-316 



2i4m 



sin 



AAm— ft sin 



4i4m — ft sin 

ft J 

7-3 1 ,must satisfy the condition: 
ft 

A x m > - 
4 

If we consider, that,d -» ; 7-31 reads: 
_ 12m 3 k 2 H 4 a 2 

P\d->0 - 7T f ) 2r T~. 7 \ AAm] 

nn u u Am _ ft sm __i 

Now we define the continuity equation: 
dP _ 

_ + vj = o 



8m 



1-n 



AAm—hsin- 



AAm 



] 7 " 



32 



7-33 



7-34 



7-35 



dP _ 32 k 2 m 2 n(n - l)t 
at [44|fc|m- ftsin 



ift 


r 3^* 






da 


2m 


r- 


l)t (n " 


i) 




4^1|fe|m" 


ftt 


cos 



_ IJJ* 



o 



da 

2An\k\mt (n - 1) 



2An\k\mt 
sin - — 



7-36 

(n-1) 



= 



Now we can solve for (t),(n).(appindixI-4); 
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In- 



±nh 



A | 4A\k\nm , 

t(n) = exp — I ; t = 

n — 1 



1.8- 
1.6- 
1.4- 
1.2- 
t 1 - 
0.8- 
0.6- 
0.4- 
0.2- 

0" 



7-37 
Live 



y 



y 



2e-07 



4e-07 



6e-07 



8e-07 



1e-06 



fig(7 — 1): shows variation of t — n, when t -> 1, we find n is indefinit. 



n{t) 



LambertW 1 In t) 



lnt 



; n = 0; n = 1 



7-38 



7-38 gives inflation at (t=l). 

Age of our universe is 10 9 y,(appendix I),n = | , from 7-37 we find: 



AAm\k\ 
nh 



= 1.059505956 x 10 6 ;m « 2.65 x 10" 29 O'-s 2 ) 



Now we can estimate radiation dominated age, with n = | ; 

W ~ 10 5 y 

Now we find the scale factor, from 7-37 and 7-38(appendix I- 6) as: 

) 7 - 39 



And d as: 
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d(n) = An 



I In 



±nh 
4A\k\nm 



n-1 



I tJi \K TO 

exp - 

I n-1 

And acceleration is given by( appendix I- 7): 

/, ±nh 

I In- 

d(n) = An(n — 1) 



exp 



I n 



4A\k\nm 



H(ri) = n 



In 



exp 



±7Tft 

4A\k\nm 

n-1 



n— 2 



appendix 1-8 to calculate mass density and Hubble constant . 
From 7-24 and( appendix 1-6) ,we can estimate the value of 

. _i 4,4m -ft sin ^ jt , 4Am-hsin^- 

sin J = ' l /_ ; or 2- = i 

\ 8m / Z 8m 



7-40 



7-41 



7-42 



7-43 



8. Complex solution: 



Now we solved 7-37 in context of the term 



nh 



4) we find: 



4A\k\m 



, according to (appendix I- 



nh 



= -4.4 x 10" 7 ± 7.6 x 10" 7 i 



8-1 



AA\k\m 

If we assume that, the constant ,A, is a complex, using 8-1, we find: 

A = (l + 1.7321i )tt, ^1* = (1 - 1.7321i)rc 8-2 
4i4m — /i sin 



8m 



= -.1276 + .2265i 



8-3 



By substituting 8-1 and 8-2 in 7-39( appendix 1-9) ,7-40( appendix I-10),7-41( 
appendix 1-1 l)and 7-42( appendix 1-12). 

From appendix 1-11, our today universe acceleration (n = |) can be assigned to the 
imaginary part of d. Also we can find the density, by substituting 8-3 in 7-31, ( 
appendix 1-13). 

If we treat the cosmic back ground radiation as black body radiation, and 
according to Stefan's law of radiation: 

oT 4 oc E 



T = 



8-4 



Where / is the ratio of redshifted radiation density. 

<f < 1 

8-3 represents maximum(limit) CBR temperature when / = 1 

If we substitute the values from which we found ,from our work, we get: 

/ = 1 
T = 5.3-3.1ifc 

When / = 0.07 

T = 2.726 - 1.58£ k 
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Now if we choose the scale factor with variable A as: 

a = A n t n 8-5 

Now we can rewrite our equations as: 

8|fc|m 7 [2\k\m \ 

P = H = sin^ — - — d 8 — 6 

4A-\k\m-hsin 4 -^ V » J 



3|fejm ^ . 2 (2\k\ma\ _ 



(ln-=^ — nln/l\ 
m ™_ J;t = 8-9 

LambertW ( ~^ t In i4t) 

n(t) = — ;n = 0;n=l 8-10 

lni4t 

T = (1.89654434 - 1.8965443381i) K 
With absolute valuej = 2.682118726 K 
When the age of universe is: 

t = 13.75 Gy 
Which is agree with recent observations. 



9. Results and Discussion: 

We have seen that the canonical quantization formalism adopted in this work 
has brought interesting results. We have solved the Wheeler-De Witt equation and 
obtained the wave function of the universe. The wave function vanishes when 
d -» . We found that density parameter is related to probability. Moreover, we have 
found that the scale factor is imaginary. From Stefan's law we have found the cosmic 
microwave background radiation (CMBR) temperature. 

Considering the model of d it is important to note that in case of a -» oo , d « 1 , i.e. 
when the universe radius is large, the expansion rate is slow, equation (2-26) which is 
the expression for the density reduces to that of GR, the density of universe decreases 
with time as radius increases in conformity with observations. It also indicates that the 
pressure, P, vanishes, equation 2.28, which also agrees with GR and observations at 
matter era. 

Thus, this case represents matter era in GR. When the radius is assumed to be very 
small with high expansion rate i.e. a -» , d -» oo, as happens at early universe, the 
pressure becomes negative , equation 2.32. This also agrees with inflation assumption. 
The model for d 2 gives also some interesting results, it shows that p oc a -3 and p = 
which can suitably describing matter era only. 

The boundary term a for this model coincides with the one emerging from that of the 
model d, as shown by equation 2.9. 
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The two cosmological constant A 1 and A 2 for the model a and d 2 show also some 

i 

interesting properties. They both vanishes at radiation era, when a = AV- . 
As seen in equations (6-2), (6-7), (6.9) and (6-13), equation (6-19) shows that the total 
cosmological constant decreases as time increases. This agrees with inflation models 
at the early universe. 

The wavefunction ipi and xp 2 of the universe for both models are also obtained in 
section [7]. The complex form of the solution in (7-13) shows the existance of 
quantum effects. The complex nature of cosmic scale factor shown in (8-2) confirms 
the existance of quantum effects. The CBR temperature is found to agree with the 
experimental range , equation 8-5. 

APPENDIX 

1 . This program (in maple code) to calculate the total cosmological constant,here 
we referred to A tot as (ty). 

>ty:=proc(t,n,k) ; 
> localresult ; 

> result : 

> 2*m*k A 2/A/ (t A n) -2*m* (A*n*t A (n-1) -t A (-n+1) /A/ (n- 

1) *k+t A (n-1) *A/ (n-1) *n A 2) * (3*A A 2*n A 2*t A (n-1) * (n-1) *t A (n- 

2) +A A 2*t A n*n* (n-1) * (n-2) *t A (n-3) ) -2*m* (A*n* (n-1) * (n- 

2) *t A (n-3) + (2*A A 3*t A n*n A 2*t A (n-1 ) * (n-1 ) *t A (n-2) -A*n*t A (n- 
1) *k-A A 3*n A 3* (t A (n-1) ) A 3) / A A 2 / (t A n) A 2) * (A A 2 * t A n*n* t A (n- 
1) +k*t) ; 

> end; 

ty :=pro((?, n, k) 

localresult, 
result; 

2xmxk A ~2/(Axt A n) - 2xmx 

(Axrixt^n - 1 ) - t\-n + 1 )xk/(Ax(n - 1 )) + t A (n - 1 )xAxn A 2l(n - 1 ))x 
(3xA A 2xn A 2xt A (n - l)x(» - 1 )xt A {n - 2) 
+ A A 2xt A nxnx(n - 1 )x(n - 2)xt A (n - 3)) - 2xmx( 
Axnx(n - l)x(n -2)xt\n -3) + ( 

2xA A 3xt A nxn A 2xt A (n - l)x(« - \)xt A (n -2) -Axnxt^n -\)xk 
-A A 3xn A 3x(t A (n - 1 )) A 3 )/(A A 2x(t A n) A 2))x(A A 2xt A nxnxt A (n -\)+kxt) 

end proc 

>ty(t,n,k) ; 
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2mk 2 



At" 



— 2 m 



Ant 



(-n + l) 
(n-l) t K 



A(n-l) 



t An 
n-l 



(3 A 2 n 2 t (n ~ 1) (n-l)t (n ~ 2) +A 2 f n (n - 1 ) (n -2) Z"~ 3) ) -2 m 



,4« (« - 1) O -2) f 



0-3) 



+ 



2^V» 2 /"" l) (»-l)/"" 2) -^ii/"" 1) Jt-^ii 3 ( "" 1) ) 



3 \ 



A 2 (t") 



(A 2 t n nt {n ~ l) +kt) 



>R0 := simplify (2*m*k A 2/A/ (t A n) -2*m* (A*n*t A (n-l) -t A (- 

n+1) /A/ (n-l) *k+t A (n-l) *A/ (n-l) *n A 2) * (3*A A 2*n A 2*t A (n- 

1) * (n-l) *t A (n-2) +A A 2*t A n*n* (n-l) * (n-2) *t A (n-3) ) - 

2*m* (A*n* (n-l) * (n-2) *t A (n-3) + (2*A A 3*t A n*n A 2*t A (n-l) * (n- 

1) *t A (n-2) -A*n*t A (n-l) *k-A A 3*n A 3* (t A (n- 

1) ) A 3) / A A 2 / (t A n) A 2) * (A A 2*t A n*n*t A (n-l) +k*t) ) ; 



RO ~-2(l0A 4 n 4 t (3n - 4) -l3A 4 n't (3n - 4) + 2A 2 n 3 t in - 2) 



k-\0A 2 n 2 t (n ~ 2) k 



A *4 2 X3«-4) (-„) 

+ 4 A n t 



r" nk 2 +4A 2 nt n ~ 2) k-t n) k 2 )mlA 



> Rl := combine (RO , power ) ; 

Rl :=-2(10A 4 n 4 t (3 "~ 4) -UA 4 ^ T" "+2A z n s r "k-\ 

+ 4A 4 n 2 P "- 4) - t ( ~ n) nk 2 +4A 2 n f ( "~ 2) k - t"' k 2 ) ml A 
> ty (t, 0, 0) ; 

> ty (t, .5,0) ; 

>ty(t, 1/2,-1) ; 




0. 



2 m 

Aft 



2 m 



3 A 2 



3 A 



8 t 



(3/2) 



A 
2Jt 



8 t 



(5/2) 



A 2 t 



I 2 



>R18 := simplify (2*m/A/t A (1/2) -2*m* (3/8*A/t A (5/2) + (- 
3/8*A A 3/t A (3/2) +l/2*A/t A (1/2) ) /A A 2/t) * (l/2*A A 2-t) ) ; 

m(6t-A 2 ) 



>R19 



expand (Rl 8 ) ; 



R18 :-- 



R19 



2t (m A 



3 m 



m A 



> ty (t, 1/2,1) ; 

2 m 



AJt 



2 m 



2t 



3 A' 



(3/2) 



A 



3 A 



St 



(5/2) 



/ 3/2) 2^ 
A 2 t 



U 2 



+ t 
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>R20 := simplify (2*m/A/t A (1/2) -2*m* (3/8*A/t A (5/2) + (- 
3/8*A A 3/t A (3/2) -l/2*A/t A (1/2) ) /A A 2/t) * (l/2*A A 2+t) ) ; 

m (6 t +A 2 ) 



R20 :-- 



2t im A 



>R21 



expand (R20) 



3m mA 

R21 := 1= + 



2f 



(3/2) 



> ty (t, 2/3,-1) 



2 m 



3 t 



2A 3 ? 

(1/3) ^ 



(1/3) \ 



,4* 



(2/3) 



(5/3) 



2 m 



27? 
+ 



8 A 



21 1 



(7/3) 



2 ?' 3/ 1/3) 

^2/4/3) 



2^ 2 / 1/3) 

t 



>R22 := simplify (2*m/A/t A (2/3) +8/27*m* (-2/3*A/t A (1/3) - 

3*t A (l/3)/A)*A A 2/t A (5/3)-2*m*(8/27*A/t A (7/3)+(- 
16/27*A A 3/t+2/3*A/t A (l/3))/A A 2/t A (4/3))*(2/3*A A 2*t A (l/3)- 
t) ) ; 



>R2 3 



R22 := 

expand (R22 ) ; 



2m(-\35t i5 ' 3) -8A 4 t im) + 96A 2 t) 



R23 :-- 



10 m 



Sit™ A 



16 mA 3 64 in A 



3 At 



(2/3) 



81 r 



27 ? 



(4/3) 



> ty (t,2/3, 1) 



8m 



2 m 



2.4 3 1 

+ 



(1/3) \ 



3 ? 



(1/3) 



At 

f 



(2/3) 



2 m 



%A 



27 

16A 3 2A ^ 



27 t 



(7/3) 



W 4/3) 



2^ 2 / 1/3) 

= + t 



J 



J 



>R24 := simplify (2*m/A/t A (2/3) +8/27*m* (- 
2/3*A/t A (l/3)+3*t A (l/3)/A)*A A 2/t A (5/3)- 
2*m* (8/27*A/t A (7/3) + ( -1 6/2 7 *A A 3 / t- 

2/3*A/t A (l/3))/A A 2/t A (4/3))*(2/3*A A 2*t A (l/3)+t)); 



2m(135* +8A*t +96AU) 



>R2 5 



expand (R2 4 ) ; 

i?25 := 



10 m 



8l/ 7/3) ^ 



16 mA~ 



3 At 



(2/3) 



81 t 



(Am A 

(4/3) 

27 t 
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>ty(t,n,0) ; 

(«-i) * v4« 
-2m Ant H , 

(3^ 2 « 2 /" _1) (n-l)/^ 2) +^ 2 ^n(n-l)(n-2)/^ 3) )-2m 
^4 « (« - 1) (n -2) t H ; 



,4 2 (0 



^ 2 f" 



>ty(t,l/2,k) ; 



2mr 



— 2m 



3* 



Ak 



3 A 8/ 3/2) 2^ 
+ 



St 



(5/2) 



^ 2 ? 



(A' 



+ kt 



>R0 := simplify (2*m*k A 2/A/t A (1/2) -2*m* (3/8*A/t A (5/2) + (- 
3/8*A A 3/t A (3/2) -l/2*A/t A (1/2) *k) /A A 2/t) * ( 1 /2 *A A 2 + k* t ) ) ; 



RO :-- 



> Rl := expand (RO) ; 



> limit (Rl, t=inf inity) ; 



k(6kt+A 2 ) 
2t m) A 



3 m k m kA 




2. This is maple code to calculate A , fi A and p A 

> lmabda=m*k A 2/A/ (t A n) -2*m* (A*n*t A (n-1) -t A (-n+1) /A/ (n- 

1) *k+t A (n-1) *A/ (n-1) *n A 2) * (3*A A 2*n A 2*t A (n-1) * (n-1) *t A (n- 

2) +A A 2*t A n*n* (n-1) * (n-2) *t A (n-3) ) +m*k A 2/A/ (t A n) - 

2*m* (A*n* (n-1) * (n-2) *t A (n-3) + (2*A A 3*t A n*n A 2*t A (n-1) * (n- 

1) *t A (n-2) -A*n*t A (n-1) *k-A A 3*n A 3* (t A (n- 

1) ) A 3) / A A 2 / (t A n) A 2) * (A A 2*t A n*n*t A (n-1) +k*t) ; 



Imabda 



2mk 2 
At" 



■ 2 m 



A n t 



(-n + 1) 
(n-1) t k 



a 2^ 

t A n 



A(n-l) n-l 



(3 A 2 n 2 t in ~ l) (n-\)t {n ~ 2) +A 2 f n (n - 1 ) (« -2) /*~ 3) ) -2 m 



(n-2) 



(»-3V 



An(n-\){n-2)t 



(n-3) 



+ 



2A'fn 2 r- L, (n-\)t {n - 2) -Ant (n - V) 

2 r +n - 



3 \ 



(n-1) 



{A 2 i*nt ~ ' + *o 



>R16 := rhs (lambda = 2*m*k A 2/A/ (t A n) -2*m* (A*n*t A (n-1) - 
t A (-n+1) /A/ (n-1) *k+t A (n-1) *A/ (n-1) *n A 2) * (3*A A 2*n A 2*t A (n- 
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1) * (n-1) *t A (n-2) +A A 2*t A n*n* (n-1) * (n-2) *t A (n-3) ) - 
2*m* (A*n* (n-1) * (n-2) *t A (n-3) + (2*A A 3*t A n*n A 2*t A (n-1) * (n- 
1) *t A (n-2) -A*n*t A (n-1) *k-A A 3*n A 3* (t A (n- 
1) ) A 3) / A A 2 / (t A n) A 2) * (A A 2*t A n*n*t A (n-1) + k*t) ) ; 

2mk 2 



R16 :-- 



At" 



— 2 m 



(n-i) t k t An 

Ant 



V 



A(n-l) n-l 



) 



(3 A 2 n 2 y~ l \n-\)t in ~ 2) +A 2 f n (n - 1 ) (n - 2) /"~ 3) ) - 2 m 

(n-3) 



An(n-\)(n-2)t 



+ 



-3 ^..2 ^+"^-Ant (n - l) k-A^n\t { "- l) ) 



2A 3 fn 2 r~ L, {n-\)t 



3 \ 



A 2 (t n ) 



(A 2 fnt (n ~ l) +kt) 
>R17 := simplif y (R16) ; 

R17:=-2(l0A 4 n 4 t i3n - 4) -l3A 4 n't <3n - 4) + 2A 2 n't (n - 2) k-l0A 2 n 2 t {n - 2) k 

+ 4A 4 n 2 P "- 4) - t n) nk 2 + 4A 2 nt ( "- 2) k- t™ k 2 ) ml A 
> R5 := expand (R17) ; 

20mA 3 n 4 (t") 26mA 3 n 3 (t") 4mAn 3 t"k 20mAn 2 t"k 

R5 := , + , v J z + 



8mA 3 n 2 (t") 2 m n k 2 8mAnt"k 2 m k 2 
t 4 + At" t 2 + At" 



> R6 := collect (R5, k) ; 



R6 ■:- 



( 2 m n 

+ 



At" A t" 



k 2 + 



4mAn 3 t" 20 m A n 2 t" 8 m Ant 



+ 



20 m A 3 n 4 (t") 26 m A 3 n 3 (t") 8 m A 3 n 2 (t") 
? + ? ? 

This is fi A 

> omega:=R17*t A 2/ (3*n A 2) ; 

a:=-2(l0A 4 n 4 t i3n - 4) -l3A 4 n 3 t (3 °- 4} + 2A 2 n 3 /"" 2) k - WA 2 n 2 k 

+ 4A 4 n 2 t {3n - 4) -t { - n) nk 2 + 4A 2 nt ( "- 2) k-t n) k 2 )mt 2 /{3An 2 ) 
> R7 := collect (omega, k) ; 
2{-t ( - n) n-t ( - n) )mt 2 k 2 



R7 :=- 



3 An' 



2(2 A 2 n 3 t { "- 2) -\0 A 2 n 2 t {n ~ 2) + 4 A 2 nt { "- 2) )mt 2 k 

3 An 2 

2(10 A 4 n 4 t (3 - 4) -13 A 4 n 3 1 (3 - 4) + 4 A 4 n 2 t (3 "- 4) )mt 2 



3 An' 



Here we referred to p A as (q): 
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> q:=R17/ (8*pi*G) ; 

q ~-(l0A 4 n 4 t i3n - 4) -l3A 4 n 3 t i3n - 4) + 2A 2 n 3 t (n - 2) k-l0A 2 n 2 t { "- 2) k 

+ 4A 4 n 2 t (3n - 4) -t { -" ) nk 2 +4A 2 nt ( "- 2) k-t ( - n) k 2 )m/(4AnG) 
> limit (R17, n=0) ; 

2k 2 m 



3. This program to calculate the wave function xp 

> wave : =proc (b, A, x) ; 
> localresult ; 

> result :psi (x) =2*sqrt (b/ (2*A*b-sin (2*A*b) ) ) *sin (b*x) ; 

> end; 

wave := p ro <(b, A, x) 

localresult, result; = 2xsqrt( b/(2xAxb - sin(2x^4x6)))xsin(6xx) 
end proc 

> wave (b, A, x) ; 

v|/(x) = 2 A / - - u — „f M< -^ A \ sin(6 x) 



2Ab-sin(2Ab) 



> wave (2*abs (k) *m/h, A, x) ; 
= 2/2 




/?7 



2 I A: I m x 

sinl 

m ^ { h 



, £ m . f 4^4 A: 

*[~ * — sin l^r 

> RO := rhs (psi (x) = 

2*2 A (1/2) * (abs (k) *m/h/ (4*A*abs (k) *m/h- 

sin (4*A*abs (k) *m/h) ) ) A (1/2) *sin (2*abs (k) *m/h*x) ) ; 

2 k mx 



R0:=2J2 




k m 



, f4^4 it m . f4A k 
h I , si] 



h 



m 



sin 



h 



> probability :=P= (RO) A 2; 

probability := P - 



8 I k I m si 



. f 2 A /» .v 



4A k 



h 



. (4A k 
- si 



m 



h 



>probability := {P = 

8*abs (k) *m*sin (2*abs (k) *m/h*x) A 2/ (4*A*abs (k) *m- 
sin (4*A*abs (k) *m/h) *h) } ; 



probability : =< P 



8 I k I m sinl 



4 A I k I m - si 



2 I A* I m x 
h 

4A\k 



m 



h 



>sl :=limit (8*abs (k) *m*sin (2*abs (k) *m/h*x) A 2/ (4*A*abs (k) *m 
-sin (4*A*abs (k) *m/h) *h) , k=0) ; 
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si :-- 



3x z 
A 3 



> limit (si, x=A) ; 



3^ 
A 



4. This program to calculate time and n from continuity equation. 

> 32*m*abs (k) *abs (k) *m*n* (n-1) * (t A (n- 

1) ) * (cos (abs (k) *m*A*n*2* (t A (n- 

1) ) /h) *sin (abs (k) *m*A*n*2* ( t A (n-1 ) ) /h) ) =0 ; 



32 m 2 | k\ 2 n (n - 1 ) t" ''cos 



2 | k | mA n t 



(n-\)\ 



sin 



2 | k | m A n t 



(n-\)\ 







>R1 := solve ( {32*m A 2*abs (k) A 2*n* (n-1) *t A (n- 

1) *cos (2*abs (k) *m*A*n*t A (n-1) /h) *sin (2*abs (k) *m*A*n*t A (n- 

1) /h) = 0}, {t}) ; 



Rl := {/ = 0}, 



^ \k\mAn ) \ \k\m An ) 



t = e 



\k\m An 
n - 1 



A: | m A n 
n-1 



>R0 := solve ( {32*m A 2*abs (k) A 2*n* (n-1) *t A (n- 

1) *cos (2*abs (k) *m*A*n*t A (n-1) /h) *sin (2*abs (k) *m*A*n*t A (n- 

1) /h) = 0}, {n}) ; 

1 \n(t)%th\ 



R0 := {n = 0}, {n = l}, 



LambertW 



4 A: 



ln(0 



1 ln(0 nth}} 
LambertW — - 



4 k mA 



ln(0 



> ty : =proc (n) ; 
> localresult ; 

> result: exp (In (c/n) / (n-1 ) ) ; 
> end; 

ty := proc(n) localresult ; result ; exp(ln(c/«)/(« - 1 )) end proc 

> ty (n) ; 



>R4 :=exp(ln(l/2*pi/ (595275 . 3 94 5*pi) /n) / (n-1) ) ; 

.-6>> 



R4 :=e 

>rl:=ty (2 /3 ) =5* 1 A 17 ; 



0.83994736660 



n-1 



rl 



8 



27c- 



=5000000000000000( 



> evalf (c) ; 
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{c = 0. 8399473661 0" 6 }, {c = -0.41997368330^ + 0.72741575760" 6 /}, 

{c^-0.41997368330" 6 - 0.7274 1 575760" 6 /} 

>r2:=ty (1/2) ; 



> tim:=l/ (4* . 83 99473 667e-6 A 2) ; 

tim := 0.3543527952 10 12 

When t= 13. 75 Gy ; 

{c = -0.44029756340" 6 -0.76261775030" 6 /}, 

{c = -0.44029756340" 6 +0.76261775030" 6 /}, {c = 0.880595 12680" 6 } 



5. This program to calculate (n) 

> np : =proc ( t ) ; 
> localresult ; 

> result : 1 /In (t) *LambertW (. 8 3 9 9473 6 6 6e-6*ln (t) *pi/pi*t) ; 

> end; 

np :=pro((0 

localresult, result; Lambert Wv0.8399473666*10 A (^KO x O/ln(O 
end proc 

>np (t) ; 

Lambert Wy0.8399473666Q- 6 ln(Q t) 
W) 

> np (10 A 12) ; 

Lambert^ 83 9947 .3 666n( 100000000000))) 
ln(100000000000p 

> evalf (%) ; 

0.5175266154 

> np (5el7) ; 

0.6666666668 

>R1 := convert (. 6720473448, ' rational ',' exact ') ; 
>R0 := convert (. 6720473448, ' rational ') ; 
>np(10 A (-36) ) ; 

Lambert\\(-0.839947366a0^ 2 rn(100000000000000000000000000000000p( 
ln( 100000000000000000000000000000000^)000 

> evalf (%) ; 

0.8399473667 10 A2 

> n=np (t) ; 

Lambert W(0. 83994736660^ ln(0 t) 

n = . . . 

ln(0 
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> np ( . 9999999999); 
> limit (np (t) , t=l) ; 



0.8399473665 10" 



0.8399473666 10" 6 

> limit (np (t) , t=inf inity) ; 

1. 

> np ( 1 ) ; 

Error, (in np) numeric exception: division by zero 

> np ( ) ; 

Error, (in In) numeric exception: division by zero 



> np (5 . 39e-44) ; 
>np ( . 9) ; 

> np (le-37) ; 
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0.4527316303 10 
0.7559526901 10" 6 



0.8399473666 10 



-43 



6. This program to calculate the scale factor (a) 

> scale: =proc (n) 
> localresult ; 

> 

> result :A*exp (In (l/2*pi/ (595275 . 3 945*pi) /n) / (n-1 ) ) A n; 

> end; 

scale :=pro<(«) 

localresult, result; ^xexp(ln(0. 8399473666* 1 A (76))/(« - 1)) A « 
end proc 

> scale (n) ; 

n 

f ( ( 0.83994736660" 6>l ^ 



A 



n 

n - 1 



> scale (1/2) ; 

> scale (2/3) ; 



595275.3973 A 



0.6299605277 10 12 ,4 

> scale ( . 66667*10 A (-6) ) ; 

0.9999998460 A 
0.18617 1 8133 1 55 y4 
0.5145804043 10^ 7 ^ 



> scale ( . 9) ; 

> scale (1.1) ; 
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This to calculate d : 

> vsc : =proc (n) ; 
> localresult ; 

> result :A*n*exp (In (1/2 *pi/ (595275 . 3 94 5*pi) /n) / (n-1 ) ) A (n- 
1) ; 

> end; 

vsc :=pro((n) 

localresult, result; ^xnxexp(ln(0.8399473666* 10 A (76))/(« -1))^ - 1) 
end proc 

> vsc (n) ; 



A n 



0.83994736660"' 



n-1 



0.8399473625 10" 6 A 
0.8399473643 10" 6 A 
0.8399473647 10" 6 A 



> vsc (1/2) ; 
> vsc (1/3) ; 
> vsc (2/3) ; 

> vsc (1/4) ; 

0.8399473685 10" 6 ^ 

> vsc (1) ; 

Error, (in vsc) numeric exception: division by zero 

> vsc (0) ; 

Error, (in vsc) numeric exception: division by zero 



> vsc ( . 66667*10 A (-7) ) ; 



> vsc ( . 99) ; 



0.8399473664 10" 6 A 



0.8399473667 10" 6 ^ 



> ad=n*exp (In (l/2*pi/ (5 95275 . 3 945*pi) /n) / (n-1 ) ) A (n-1 ) ; 



ad = n 



0.83994736660 



-6 



n-1 



Y 



> vsc ( . 98) ; 



> limit (vsc(n) ,n=l) ; 



> limit ( vsc (n) , n=0 ) ; 



0.8399473654 10" 6 ^ 



0.8399473666 10" 6 ^ 



0.8399473666 10" 6 A 



> limit ( vsc (n) , n=0 ) ; 



0.8399473666 10" 6 ^ 
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7. This program to calculate the acceleration 

> acc : =proc (n) ; 
> localresult ; 
> 

> result :A*n* (n-1) *exp (In (l/2*pi/ (595275 . 3 94 5*pi) /n) / (n- 
1) ) A (n-2) ; 

> end; 

acc := pro ((«) 

result; 

Axnx(n - 1 )xexp(ln(0. 8399473666* 10 A (76))/(« - 1 )) A (« -2) 
end proc 

> acc (n) ; 



An(n - 1 ) 



0.83994736660 



n - 1 



> acc (1/2) ; 

> acc (2/3) ; 

> acc ( . 99) ; 

> acc (le-36) ; 



-0.1185185168 IQ- 17 A 



-0.5599649062 10" 24 ^ 



-0.6108574821 10^ 15 ^ 



-0.7055115763 10 24 A 



> acc (1) ; 

Error, (in asc) numeric exception: division by zero 
> acc (0) ; 

Error, (in asc) numeric exception: division by zero 



> acc (-1 ) ; 

> acc (1 . 1) ; 



(0.9473319314 10 18 -0.1539600709 10 * I) A 



0.1246402218 10 55 ^ 

> acc ( . 83 99473 666e-6) ; 

-0.8399466611 10" 6 A 

8. This program to calculate density: 

> rho : =proc (n) ; 
> localresult , • 

>result:3.61e9* (n/ (exp (In (l/2*pi/ (595275 . 3 945*pi) /n) / (n- 
1) ) ) ) A 2* (sin (1 . 8 7e6*n*exp (In (l/2*pi/ (595275 . 3 945*pi) /n) / ( 
n-1) ) A (n-1) ) ) A 2/A; 

> end; 
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p :=pro((«) 

localresult, 
result; 

O.361*lO A 10« A 2x 

sin(0.187*10 A ^jxexp(ln(0.8399473666*10 A (76))/(n - \))\n - 1)) A 2/ 
(exp(ln(0.8399473666*10 A (y6))/(« - l)) A 2x^) 

end proc 

> rho (n) ; 



0.361 10 1( V sin 



0.187 10 7 n 



0.83994736660' 



n - 1 



■6 \\ 



> rho (2/3) 



0.83994736660"' 



n - 1 



A 



0.64177776340 



■26 



> rho (1/2) 



> rho ( . 9) ; 



0.7187460523Q- 14 
A 



0.73482426110 



i i i 



>rho(.8399473668e-6) ; 



0.00254689677' 



>rho(.4527316303e-49) ; 



0.00254689678' 



Hubble parameter: 

> localresult ; 
> 

> result :n/ (exp (In (l/2*pi/ (5 95275 . 3 945*pi) /n) / (n-1) ) ) ; 

> end; 

hu :=proc(n) localresult ; resw/f ; n/exp(ln(0. 8399473666* 10 A (-6) /n)/(n - 1 )) end proc 

> hu (n) ; 



0.83994736660 



« - 1 



■6 A A 



>hu (1/2) ; 



0.1411023144 10 



■i i 
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>hu (2/3) ; 



> hh:=H=hu (n) ; 



0.1333333325 10" 17 



hh :=H = 



0.83994736660 



■6 ^ 



n-1 



> limit (hh, n=0) ; 



> limit (hh, n=l) ; 



> hu ( . 9 ) ; 



>hu( . 001) ; 



H= 0.8399473666 10" 



i/ = Float( undefined ) 



0.4511678498 10 



-60 



0.6618041261 10" 



redshift: 

> z : =proc (nO , n) ; 
> localresult ; 

> 

> result: (exp (In ( ( . 83 99473 667e-6) /nO) / (nO- 

1) ) ) A nO/ (exp (In ( ( . 83 99473 667e-6) /n) / (n-1 ) ) ) A n; 

> end; 

z :=pro<(«ft, n) 

localresult, 
result; 

exp(ln(0.8399473667 5|: 10 A (V6P)/(n0-l)) A n0/ end proc 
exp( ln(0. 8399473667* 10 A (V6))/(« - 1 )) A n 



> z (nO, n] 



^ ( fo.83994736670" 6 % 
H nO 
nO-i 



nO 



> z (2/3, 1/2) ; 
> z (2/3, 1/3) ; 



0.83994736610 



■6 ^ 



n-1 



0.1058267368 10 7 



0.1000000003 10 



10 



> z (2/3, . 83 99473 6 66e-6) ; 

0.6299605277 10 12 



>z(.69 62 745918,.839 94 7366 6e-6); 

0.3697571388 10 



14 
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9. Imaginary scale factor: 

> scale : =proc (n) ; 

> localresult; 

> 

> result: (-3.13999999 9-5.438 63 9534*1)* (exp(ln( (- 
.44 02 975 634e-6+.7 62 6177503e-6*I) /n) / (n-1 ) ) ) A n; 

> end; 

scale :=proc(«) 
localresult; 
result; 

( - 3.139999999- 5.438639534*1) * (exp(ln( ( - 4.402975634 
•10 A -7 + 7.626177503 10 A -7*I) ln)l(n - l))) A n 
end proc 

> scale (n) ; 

(-3.139999999 

-7, ^ \ n 



- 5.4386395341) 

> scale (1/2) ; 

> scale (2 /3 ) ; 



In 



-4.402975634 10"' + 7.626177503 10"' I 



1.78288518510 6 - 3.08804772110 6 I 



1.79967703510 12 - 3.1 1713206710 12 I 



> (scale (2/3) /scale (1/2) ) A 3*4/3*3 . 14159/5 . 67e-8 ; 

-7.598341041 10 25 - 5.22670213810 16 I 



10. This program to calculate imaginary d : 

> vsc : =proc (n) ; > vsc : =proc (n) ; 

> localresult; 

> 

> result: (-3.13999999 9-5.438 63 9534*1) *n*exp(ln( (- 
.44 05 62 6353e-6+.7 6307 68 683e-6*I) /n) / (n-1 ) ) A (n-1 ) ; 

> end; 

vsc :=proc(«) 

localresult; 
result; 

( - 3.139999999- 5.438639534*1) *n* (exp(ln( ( - 4.40562635 
■10 A -7 + 7.630768683 10 A -7* I) ln)l(n - 1))) A (« - 1) 
end proc 

> vsc (n) ; 
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(-3.139999999 

- 5.4386395341) 

, n — 1 



In 



-4.405626353 10"' + 7.630768683 10"' I 



-0.000005533466694+ 3.01816685610" 15 I 



> vsc (1/2) ; 
> 

> vsc (2/3) ; 

> 

-0.000002766733355- 0.000004792122748 

> vsc ( . 8399473666e-6) ; 

0.000005533466698- 2.95657499510" 15 I 

> vsc ( .4527316303e-49) ; 

0.000005533466678- 2.32253429910" 15 I 

> vsc ( . 9) ; 

-0.000001709935257- 0.000005262639581 

> limit (vsc (n) , n=0 ) ; 

0.000005533466698- 1.531 19181510" 16 I 

> limit (vsc (n) , n=l) ; 

-0.000002766733348- 0.00000479212273Q 



1 1 . This program to calculate imaginary a 

> ac : =proc (n) ; 

> localresult; 

> 

> result: (-3.1415 9-5.4413 93498*1) *n* (n-1) *exp(ln( (■ 
.44 05 62 6353e-6+.7 6307 68 683e-6*I) /n) / (n-1) ) A (n-2) ; 

> end; 

ac ~proc(n) 

localresult ; 
result; 

( - 3.14159- 5.441393498*1) *n*(n - l)*(exp(ln(( 
- 4.405626353 10 A -7 + 7.630768683 10 A -7*I) ln)l{n - 1))) 
A (« - 2) 
end proc 



> ac (n, 



(-3.14159- 5.4413934981) n (n 



-4.405626353 10" 7 + 7.630768683 10" 7 I 



x«-2 



n - 1 



> ac (1/2) ; 



-4.29825792010" 



7.44480108210" 18 I 
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> ac (2/3) ; 



2.13035835110" 24 + 3.68988 892410" 24 I 



> ac ( . 8399473666e-6) ; 

0.000002903846946- 0.000005029590028 

> 

> ac ( .4527316303e-49) ; 

5.38745900010 37 - 9.33135273810 37 I 

> limit (ac (n) , n=0 ) ; 

Float (undefined ) 

> limit (ac (n) , n=l ) ; 

Float (undefined ) 



12. This program to calculate imaginary H 

> hu : =proc (n) ; 
> localresult ; 

> result :n/exp (In ( (-. 44 05 62 6353e-6+ . 7 63 7 68 68 3e- 
6*1) /n) / (n-1) ) ; 

> end; 

hu :=pro in) 

localresult, 
result; 

n/exp(ln(-0.4405626353 *10 A (-&)0. 7630768683* \0 A (^n)/(n 
end proc 

>hu (n) ; 

n 

' f-0.44056263530" 6 + 0.76307686830" 6 /^^ 



>hu (1/2) ; 
>hu (2/3) ; 
> hu ( . 9 ) ; 
>hu ( . 8399473666e-6) ; 



0.7763817423 10" 12 -0.1344732621 10" 11 / 



0.1539201547 10" 17 + 0.2801983311 10" 26 / 



0.3640524881 10" 60 + 0.6305574119 10" 60 / 



-0.4405639950 10" 6 + 0.7630761242 10" 6 / 

>hu (. 71619913179) ; 

0.4899236296 10 21 +0.9545864800 10" 21 / 

>hu (1/3) ; 

-0.1432574160 10" 8 + 0.8449225532 10" 18 / 
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13. This program to calculate imaginary density 

> rho : =proc (n) ; 
> localresult ; 

> 

> result: (3/8 / (-.127608229 6+. 21 64 9 92 92 8*1) /3 . 14/6 . 672e- 
11) * (n/exp (In ( (-. 4 4 02 97 5 634e-6+ . 7 62 6177503e-6*I) /n) / (n- 
l))) A 2*(sin((3.14/2/ (-. 4402 975634e-6+ . 7 62 6177503e- 
6*1) *n* (exp (In ( (- . 44 02 975 634e- 6+ . 7 62 6177503e- 6* I ) /n) / (n- 
1) ) ) A (n-1) ) ) A 2) ; 

> end; 

p := pro ((h) 

localresult, 
result; 

(-0.36 16687599* 1 A H)0. 6 136048669* 10 A iC)xn A 2xsin( 
(-89 1442. 5893- 0.1 544023857* 10A7)xnx 

exp (ln( -0.4402975634* 10 A (-6) + 0.76261 77503 *10 A (-6) xl/n)/(n - 1 )) A 
(n-l)) A 2/ 

exp (ln( -0.4402975634* 10 A (-6) + 0.76261 77503 *10 A (-6) xl/n)/(n - 1 )) A 2 

end proc 

> rho (n) ; 



(-0.36166875990 10 - 0.6 136048669 10 /)« 2 sin 



/ i 



(891442.5893h0.15440238570 7 /)n 

. 2 



-0.4402975634 0" 6 +0.76261775030" 6 



n-1 



r [ -0.4402975634 0" 6 + 0.7626 1 77503 0" 6 1^ ^ 



n-1 



> rho (1/2) ; 



0.1713113633 10 - 1 3 -0.1542238948 10" 15 / 



> rho (2/3) ; 



0.5045186728 10 25 - 0.3943643036 10" 25 / 



14. Imaginary redshift: 

> z : =proc (nO , n) ; 
> localresult ; 
> 

> result: (n/nO) * (exp (In ( (- . 4 4 02 97 5 634e-6+ . 7 62 6177503e- 
6*1) /n0) / (n0-l) ) ) A n0/ (exp (In ( ( — . 4402975634e — 
6+.7 62 6177503e-6*I) /n) / (n-1) ) ) A n; 
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> end; 

z :=pro <(nO, n) 
localresult, 
result; 

«xexp(ln(-0.4402975634* 10 A (-6)0. 7626177503 * \0 A (-^n0)/(n0 - \)) A n0 
/(nOx 

exp(ln(-0.4402975634*10 A (-6p.7626177503*10 A (-^«)/(«- \)) A n) 
end proc 

> z (n0,n) ; 



f ( ( -0.4402975634 0" 6 + 0.7626 1 77503 0" 6 7 ^ 
nO 

nO- 1 



nO 



nO 



> z (2/3, 1/2) 
> z (2/3, 1/3) 



-0.44029756340" 6 + 0.7626 1 77503 0" 6 



;j - 1 



378531.8851 -655636.4568 / 



0.2326813813 10 9 - 0.4030159730 10 9 / 



> z (2/3, . 8399473666e-6) ; 

721249.4858 - 1.267934611 / 

> z ( . 69 62 745918, .839 94 7366 6e- 6); 

0.3639034496 10 8 + 0. 1653402114 10 8 / 

> z (1/2, . 8399473666e-6) ; 

280831.8725 +486417.0481 / 



> limit (z (nO, n) , nO=l) ; 
> limit (z (nO, n) , nO=0) ; 



0. 

1 



-0.44510483830" 6 + 0.770944 1 941 0" 6 



n-1. 
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